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In this work, we have considered interacting dynamical model taking variable modified Chaplygin 
gas which plays as dark energy coupled to cold dark matter in the flat FRW universe. Since the 
nature of dark energy and dark matter is still unknown, it is possible to have interaction between 
them and we choose the interaction term in phenomenologically. We have converted all the equations 
in the dynamical system of equations by considering the dimensionless parameters and seen the 
evolution of the corresponding autonomous system. The feasible critical point has been found and 
for the stability of the dynamical system about the critical point, we linearize the governing equation 
around the critical point. We found that there exists a stable scaling (attractor) solution at late 
times of the Universe and found some physical range of n and the interaction parameter c. We 
have shown that for our calculated physical range of the parameters, the Universe explores upto 
quintessence stage. The deceleration parameter, statefinder parameters, Hubble parameter and the 
scale factor have been calculated around the critical point. Finally some consequences around the 
critical point i.e., the distance measurement of the Universe like lookback time, luminosity distance, 
proper distance, angular diameter distance, comoving volume, distance modulus and probability of 
intersecting objects have been analyzed before and after the present age of the Universe. 

PACS numbers: 98.80.Cq, 98.80.-k, 95.35. +d, 95.36.+X 

I. INTRODUCTION 

Recent observations of redshift and luminosity of type la supernovae [U-IH , WMAP [5] , Chandra X-ray 
Observatory [6| etc. indicate that universe is spatially flat and undergoing accelerated expansion. The most 
important discovery over the last few decays is to search for the existence of dark energy which violates 
the strong energy condition i.e., p + 3p < [7-9]. The mystifying fluid namely dark energy is understood 
to dominate the 70% of the Universe and have enough negative pressure to drive current acceleration 30% 
dark matter (cold dark matters plus baryons). There are various candidates to play the role of the dark 
energy. The most obvious candidate for dark energy is the cosmological constant with the equation of 
state p/p = —1. The type of dark energy represented by a scalar field is often called quintessence. Other 
dark energy candidates are namely tachyonic field [loj . DBI-essence [ll|, Chaplygin gas flp . phantom, 
holographic dark energy hessence dark energy [l3| , k-essenece and dilaton dark energy [l5[ . 

Recently the Chaplygin gas, also named quartessence, characterized by an exotic equation of state 
p = —B/p (B > 0), was suggested as a candidate of a unified model of dark energy [T3 |. The Chaplygin gas 
behaves as pressureless fluid for small values of the scale factor and for large values of the scale factor as a 
cosmological constant, tends to accelerate the expansion. The above equation was generalized to the form 
p = -B/p a , < a < 1 [Hill]. Consequently this was modified top = Ap-B/p a , < a < 1, A>0 which 
is known as modified Chaplygin gas [19h21| | shows a radiation era (A = 1/3) when the scale factor a(t) is 
vanishingly small and a ACDM model when the scale factor a(t) is infinitely large. Variable Chaplygin gas 
first proposed by Guo and Zhang [22| with equation of state p = —B/p, where B is a positive function of the 
cosmological scale factor 'a' i.e., B = B{a). This assumption is reasonable since B(a) is related to the scalar 
potential if we take the Chaplygin gas as a Born-Infeld scalar field [23| • Afterward there are some works 
on variable Chaplygin gas model [2J, [2||. Further, Debnath [26| introduced variable modified Chaplygin 
gas (VMCG) for acceleration of the universe. Several interesting features and physical interpretations of 
VMCG have been shown by several authors [27M30T | . The dynamical system analysis of pure and generalized 
Chaplygin gas model and the nature of critical points have been analyzed for Einstein's gravity and loop 
quantum gravity [3ll - [36| . 
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In this work, we consider a model of interacting variable modified Chaplygin gas (VMCG) with dark 
matter in the framework of Einstein gravity. We construct the formalism of autonomous dynamical system 
of equations for this interacting model. Since the nature of dark energy and dark matter is still unknown, it 
is possible to have interaction between them and we choose the interaction term in phenomenologically. We 
convert them to dimensionless form and perform stability analysis and solve them numerically. We obtain 
a stable scaling solution (which is also an 'attractor') of the equations in FRW model. Some consequences 
like lookback time, proper distance, luminosity distance, angular diameter distance, comoving volume, dis- 
tance modulus and probability of intersecting objects of the solution around the critical point have been 
investigated. We discuss our results in the final section. 

II. DYNAMICAL MODEL OF INTERACTING VMCG AND DARK MATTER 

We consider a spatially flat universe with VMCG as dark energy and dark matter interacting through an 
interaction term. Thus the Einstein equations and continuity equation of VMCG and dark matter can be 
written respectively as (choosing 8ttG = c = 1) 

«' = i, (i) 

H = -\{P + P) (2) 

and 

p + 3H(p + p) =0 (3) 

where p = p vmcg + Pdm an d p — Pvmcg + Pdm (Pdm is very small quantity, somewhere dark matter assumed 
as pressureless quantity) are the total cosmic energy density and pressure respectively with the subscripts 
vmcg and dm denote the VMCG and dark matter respectively. Since we consider that the VMCG and dark 
matter do not conserve separately, so their continuity equations are 

Pvmcg ~t~ 3-H (pvmcg ~t~ Pvmcg) Q (4) 

and 

Pdm + 3i? (pdm + Pdm) = Q (5) 

To obtain a suitable evolution of the Universe an interaction is often assumed such that the decay rate 
should be proportional to the present value of the Hubble parameter for good fit to the expansion history 
of the Universe as determined by the Supernovae and CMB data [13] ■ Here Q is the interaction term 
which have dimension of density multiplied by Hubble parameter. For suitable choice of Q is Q = 3cHp, 
c is the coupling parameter denoting the transfer strength. The interaction term can not be trace out 
from the first principles due to unidentified scenario of both the dark energy and dark matter. If Q < 
then energy density of dark energy become negative at sufficiently early times, therefore the second law 
of thermodynamics can be violated [381 ]. Thus Q must be positive and small. Also the observational 
data of 182 Gold type la supernova samples, CMB data from the three year WMAP survey and the 
baryonic acoustic oscillations from the Sloan Digital Sky Survey estimated that the coupling parameter 
between dark matter and dark energy must be a small positive value (of the order unity), which sat- 
isfies the requirement for solving the cosmic coincidence problem and the second law of thermodynamics j39l | . 



The VMCG equation of state is given by \2 



Pvmcg — Ap vmcg with A, a are constants and < a < 1 (6) 

Pvmcg 
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where B(a) is a positive function of the cosmological scale factor 'a'. Now consider for simplicity, 
B(a) = Bga~ n , where n and Bo are constants. More precisely the restriction for accelerating universe 
becomes < n < 4 [26j . 

The dark matter equation of state is 

Pdm = WdmPdm, with w dm is a small constant. (7) 

To analyze the dynamical system, we convert the physical parameter into dimensionless form as follows: 

r — In n II — O - Pvmc 3 anr j „ — Pvmcg ,„v 

x — ma, u — siomcj — ana u — 3i/ 2 

The equation of state of the VMCG can be expressed as 

/ \ Pvmcg V /„N 

Pvmcg W 

Defining dimensionless density parameters of dark matter Q dm = ffp- and using the Friedmann equation 
(1), we obtain 

^dm — 1 ^vrncg — 1 ^ (^^) 

Thus for the flat universe, u must lies in the region < u < 1, since the energies of VMCG and dark 
matter can not be negative. 

Now we can cast the evolution equations in the following autonomous system of u and v in the form: 

du 

— = — 3c — 3 (u - 1) {uwd m - v) , (11) 



u 



*_{« 1 + ,- a £}{-*-,( 1 + £).} + „(^S) 

+3v{(l + ^)u + (l + w (im )(l-u)} (12) 
A. Critical point 

The critical points of the above system are the solution of the equations — — 0. The only feasible 
critical point is obtained as 

_ n - 3 (1 - c + w dm ) (1 + a) 

u crit — j r77" ^ ) 

n - 3 (1 + w drn ) (1 + a) 
n 2 - 3n (2 + w dm ) (1 + a) + 9 (1 + w dm + cw dm ) (1 + a) 2 
Wc " 4 ' 3(l + a){3(l + uw)(l + a)-n} [ ' 

Since in a spatially flat universe, the physically meaningful range of u is < u < 1, hence < u cr n < 1 
leads to the condition < c < 3 ( 1 + I "dmKi+a)"" a ^ on g n < m i n {3(l + w dm )(l + a), 4}, [26| for existence 
of the critical point. This condition of c represents there is an energy transfer from to dark matter to VMCG. 



B. Stability around critical point 

For the stability of the dynamical system about the critical point, we linearize the governing equation 
around the critical point i.e., u — u cr n + Su and v — v cr it + Sv, we obtain 
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du 

S I -T- ) = [3 (v + w dm - 2uw dm )] crit Su + [3 (-1 + u)] crit 5v, 



3v(c + v)a ., „,„ ., 
-3w dm v v = ' + A{n - 3(1 + a)} 



5u 



{3 — 3A — n + 6v + 3(1 — u)w dm } + 



crit 

3(c + u - Au + 2v)a 
u 



Sv 



(15) 



(16) 



We find the eigen values of the Jacobian matrix at the critical point (for simplicity we set a = 1) as in 
the following: 



3n - 6(3 + 2A + c + w dm ) + 



36c 2 



6c + n - 6(1 + w dm ) 



± 



1 /{(6 - n) 2 + 72cw d . 

iv 



36w L + 12A (6 - 6c - n + 6w dm )} 



{6c + n - 6(1 + w dm )} 



' {-144c 2 + (6 - n) 2 + 12A (6 - 6c - n + 6w dm ) - 36w 2 m - 24c (n - 6 - 9w dm )} 



{6c + n — 6(1 + w dm )} 



(17) 
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FIG. 1: The dimensionless parameters are plotted against e-folding time for different values of n. The initial condition 
is v(0) — —0.16, u(0) = 0.6. Other parameters are fixed at c = 0.01, A — 0.3, Wdm = 0.01 and a = 0.8. 

FIG. 2: The evolution of w(x) = v/u for the interacting VMCG model, corresponding to the initial conditions 
«(0) = 0.6, v(0) — —0.016; The figure is drawn for different values of n, other parameters are fixed at c = 0.01, A = 0.3, 
Wdm = 0.01 and a = 0.8 respectively. 



If the real parts of the above eigenvalues are negative, the critical point is stable node and is a stationary 
attractor solution; otherwise unstable and thus oscillatory. Physical meaningful range of c is < c < 
3(i+wdrnKt+a)-n j n ran g e t ne critical point (u cr i t ,v crit ) is stable and is a late-time stationary 
attractor solution. Here we plot some figures to show the properties of evolution of the Universe controlled 
by the dynamical system (11) and (12). The dimensionless parameters u and v have been drawn in figure 1 
in terms of x = In a. Also the EOS parameter w vmcg = ^ is drawn in figure 2 for n — 1, 2, 3. We see that u 



■5 



becomes positive but v and - explore to the negative level above —1. The phase-space diagrams have been 
shown in figures 3 and 4 for n = 1, 3 respectively. We see that from the progressions of the phase-portrait, u 
goes to 1 and v tends to a negative value above —1 and hence the solution is stationary attractor and thee 
corresponding critical point is a stable node. 
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FIG. 3: The phase space diagram of parameters depicting an attractor solution. The initial conditions chosen are 
v(0) = -0.16, u(0) = 0.6 (green); w(0) = -0.12, u(0) = 0.6 (blue); v(0) = -0.08, it(0) = 0.7 (red); v(0) = -0.4, 
w(0) = 0.8 (brown). Other parameters are fixed at c = 0.01, A = 0.3, Wdm = 0.01, a = 0.8, n = 1. 



FIG. 4: The phase space diagram of parameters depicting an attractor solution. The initial conditions chosen are 
v(0) = -0.16, u(0) = 0.6 (green); v(0) = -0.12, u(0) = 0.6 (blue); v(0) = -0.08, u(0) = 0.7 (red); v(0) = -0.4, 
tt(0) = 0.8 (brown). Other parameters are fixed at c = 0.01, A — 0.3, Wdm = 0.01, a = 0.8, n — 3. 



At the critical point, the deceleration parameter q = — 1 — yH/H 2 J is obtained as 

3 

q=-l+-X, X = 1 + V crit + W dm (l - U crit ) (18) 

Integrating, we have the Hubble parameter as (ignoring the integrating constant) 

" - ir t < 19 > 

Again integrating (19), we get 



( t \ 3X 

a(t) = a (20) 

which gives power law expansion of the universe, do is the present value of the scale factor when t = to 
is the present time which is given by to = 3XH , Since we have obtained the deceleration parameter q, 
Hubble parameter H and scale factor a(t) around the critical point, so these are valid for the late stage 
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FIG. 5: The deceleration parameter q is plotted against the evolution of the universe for different values of n. Other 
parameters are fixed at c = O.Of, A = 0.3, Wd m = 0.01, and a = 0.8. 

FIG. 6: The deceleration parameter q is plotted against the evolution of the universe for n = 0. Other parameters 
are fixed at c = 0.01, A = 0.3, Wdm = 0.01, and a — 0.8. 

of the evolution of the universe i.e., when time t is very large. Thus X should be taken positive value for 
positivity of the Hubble parameter H and hence the scale factor a(t) is increasing function of time t. So from 
(18), we obtain q G (—1, 0) for accelerating phase of the universe if X 6 (0,2/3). So the VMCG couples with 
dark matter filled in the universe always drives acceleration. Figure 5 shows the graph of q against x = In a 
for some particular values of n — 1, 2, 3. For n = 0, the VMCG model reduces to MCG and these couples 
with dark matter filled in the universe always drives acceleration and this model explores upto ACDM stage 
(q = -1) (figure 6). 

C. Statefinder parameters 

We also determine the dimensionle ss p air of cosmological diagnostic pair {r, s} dubbed as statefinder 
parameters introduced by Sahni et al [21] . The two parameters have a great geometrical significance since 
they are derived from the cosmic scale factor alone, though one can rewrite them in terms of the parameters 
of dark energy and matter. Furthermore, the pair characterize the properties of dark energy in a model- 
independent manner i.e. independent on the theory of gravity. Also this pair generalizes the well known 
geometrical parameters like the Hubble parameter and the deceleration parameter. The parameter r forms 
the next step in the hierarchy of geometrical cosmological parameters H and q. The parameters are given 
by 

-zM 1 -x) (1 - 3X) - '-=m^wr x (21) 



where X is obtained in equation (18). Since we have already obtained that for accelerating universe, 
X < 2/3. In this range, s is always positive. But r>0forA"<l/3 and r<0forl/3<X< 2/3. 
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FIG. 7: The statefinder parameter r is plotted against the evolution of the universe for different values of n. Other 
parameters are fixed at c = 0.01, A = 0.3, Wd m = 0.01, and a = 0.8. 

FIG. 8: The statefinder parameter s is plotted against the evolution of the universe for different values of n. Other 
parameters are fixed at c = 0.01, A = 0.3, Wdm = 0.01, and a = 0.8. 

III. CONSEQUENCES: DISTANCE MEASUREMENT OF THE UNIVERSE 

In cosmography (the measurement of the Universe) there are many ways to specify the distance between 
two points, because in the expanding and accelerating Universe, the distances between comoving objects are 
constantly changing, and Earth-bound observers look back in time as they look out in distance. The unifying 
aspect is that all distance measures somehow measure the separation between events on radial null trajecto- 
ries, i.e., trajectories of photons which terminate at the observer. Here we will compute various cosmological 
distance measures. In this section, we shall discuss the lookback time, luminosity distance, proper distance, 
angular diameter distance, comoving volume, distance modulus and probability of intersecting objects. 

A. Lookback Time 

The lookback time to an object is the difference between the age of the Universe now (at observation) and 
the age of the Universe at the time the photons were emitted (according to the object). As light travels with 
finite speed, it takes time for it to cover the distance related to the redshift it encountered. So, a look into 
space is always a look back in time. It is used to predict properties of high redshift objects with evolutionary 
models, such as passive stellar evolution for galaxies. Thus if a photon emitted by a source at the instant t 
and received at the time to then the photon travel time or the lookback time to — t is defined by |40H44| 



t-to 



a da 

„ a 



(22) 



where ao is the present value of the scale factor of the universe and can be obtained from (20) at t = to. 
The redshift is an important observable as they can be measured easily from the spectral lines and the 
redshift increases of an object with its distance from us. Lookback time is used to predict properties of 
high-redshift objects with evolutionary models, such as passive stellar evolution for galaxies. The redshift z 
can be defined by 



ao 
a 



l + z 



(23) 



which gives the lookback time in the following form 
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*-*o = ^ F tt nr-i ( 24 ) 



3XH, 



o 



(1 + *) 



For accelerating universe we have already get X < |. Early universe is represented by z — >• oo implies 
t — >• and late universe z — > — 1, which equivalently implied t — > oo. Also z —> gives the present age t — > tg 
of the universe. 



B. Proper Distance 



As light needs time to get from an object to the observer, one can define a distance that may be measured 
between the observer and the object with a ruler at the time the light was emitted, the proper distance. 
When a photon emitted by a source and received by an observer at time to then the proper distance between 
them is defined by HJfil H H| 

[ ao da f to dt 

d = a — = a a / — (25) 
Ja aa Jt a 

which gives 

'-ssdM 1 "^)*^} (26) 

The proper distance may also called the comoving distance (line of sight) of the Universe in today. So 
between two nearby objects in the Universe, the distance between them which remains constant with epoch if 
the two objects are moving with the Hubble flow. In other words, it is the distance between them which would 
be measured with rulers at the time they are being observed divided by the ratio of the scale factor of the 
Universe then to now. Next thing to be defined is the transverse comoving distance, which is a quantity used 
to get the comoving distance perpendicular to the line of sight. For flat universe, the transverse comoving 
distance is always identical to the comoving distance (line of sight). That means the transverse comoving 
distance = proper distance = d, for our model of the flat FRW Universe. 
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FIG. 9: Look back time, proper distance, luminosity distance and angular diameter distance as a function of redshift 
assuming Ho ~ 72km/s/Mpc before the present age of the universe. 



FIG. 10: Look back time, proper distance, luminosity distance and angular diameter distance as a function of redshift 
assuming Ho = 72km/s/Mpc after the present age of the universe. 
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C. Luminosity Distance 

If L be the total energy emitted by the source per unit time and £ be the apparent luminosity of the object 
then the luminosity distance is defined by [40T - |42l EBI . l4rj ] 

/LV .,. 2 ( ^ 1 



d '=UsJ -^-mm^X^-^W! (27) 

D. Angular Diameter Distance 

The angular diameter of a light source of proper distance D observed at to is defined by [4fj| - l43l . liH [4(| 

dkT^ (28) 

The angular diameter distance (d A ) is defined as the ratio of the source diameter to its angular diameter 
(in radians) as 

d A = - = d L {l + z)- 2 = d(l + z)- 1 (29) 
For our model angular diameter distance (d A ) is given by, 



dA -H (3X-2)\—z (l + z)^- 2 ) (30) 

It is used to convert angular separations in telescope images into proper separations at the source. It is 
famous for not increasing indefinitely as z — > oo; it turns over at z ~ 1 and thereafter more distant objects 
actually appear larger in angular size. The angular diameter distance is maximum at 



3X - 2 

and corresponding maximum angular diameter d A \ max taking the form 



O \ 3(2-X) 

" -1 (31) 



d A \ ma. 



H (3X - 2) 



3X-4J \3X-4 



2 \ 2- 



1 - I 1 \ 3 ( X " 2 ) 1 / 1 \ 3(X-2) 

2 1+ 3<x-2) — 2 < 4 3 < 2 - x ) 



(32) 



Look back time, proper distance, luminosity distance and angular diameter distance before and after the 
present age of the universe are respectively drawn in figures 9 and 10. 

E. Comoving Volume 

The comoving volume Vc is the volume measure in which number densities of non-evolving objects locked 
into Hubble flow are constant with redshift as [H, |H, liEj 

dV c = D H ( 1+ J^ dA dndz = ^-{l + zf-^d 2 A dVLdz (33) 
E(z) Ho 

where dQ is the solid angle element and d A is the angular diameter, E(z) = - and Dh = is the 
Hubble distance (c is the velocity of light) and in our model we assume c = 1, dQ = 1, H Q = 72km/s/Mpc. 
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FIG. 11: Comoving volume element dVc/dz as a function of redshift assuming Ho = 72km/s/Mpc, dQ = 1, before 
the present age of the universe. 



FIG. 12: Comoving volume element dVc/dz as a function of redshift assuming Ho = 72km/s/Mpc, dQ = 1, after 
the present age of the universe. 



So, the comoving volume is proper volume times the ratio of scale factors now to then to the third power. 
Comoving volume element dVc/dz before and after the present age of the universe are respectively drawn 
in figures 11 and 12. We see that when z > 0, the volume element decreases as z decreases and when z < 0, 
the volume element increases as z decreases. 
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FIG. 13: Distance modulus Dm as a function of redshift before the present age of the universe. 
FIG. 14: Distance modulus Dm as a function of redshift after the present age of the universe. 



F. Distance Modulus 



The distance modulus is define by 



D M = 5 log 



lOpc 



(34) 
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because it is the magnitude difference between an objects observed bolometric (i.e., integrated over all 
frequencies) flux and what it would be if it were at 10 pc (this was once thought to be the distance to Vega) 
and dh is the luminosity distance. Distance modulus Dm as a function of redshift before and after the 
present age of the universe have been shown in figures 13 and 14. We see that for z > 0, Dm decreases as z 
decreases but for z < 0, Dm increases as z decreases upto a certain stage (about z ~ —0.6) and after that 
Dm decreases as z decreases. 

G. Probability of intersecting objects 

The incremental probability dP that a line of sight will intersect one of the objects in redshift interval dz 
at redshift z is given by [HI, |43[ 

(1 + z) 2 

dP = n(z)<T{z)D H K ' dz (35) 

where n[z) is the comoving number density and tr(z) areal cross-section. Assuming n(z)a(z) = 1, we 
obtain 

dP= -^-{l + zf-^-dz (36) 
For our model the expression of Probability of intersecting objects becomes 

p -1^53Ej{( 1+ "- ¥ - 1 } (37) 

In figures 15 and 16 we draw intersection probability P as a function of redshift before and after the 
present age of the universe respectively. We see that for z > 0, P decreases as z decreases and for z < 0, P 
increases as z decreases. 
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FIG. 15: Intersection probability P as a function of redshift before the present age of the universe. 
FIG. 16: Intersection probability P as a function of redshift after the present age of the universe. 



IV. DISCUSSIONS 



In this work, we have considered flat FRW model of the universe in Cosmology. We have assumed that 
the dark energy can be considered in the form of variable modified Chaplygin gas (VMCG). The interaction 



12 



between dark matter and VMCG has been investigated in this model. To analyze the dynamical system, 
we have converted the physical parameters into dimensionless form. We have made a comprehensive 
phase-plane analysis for the VMCG model. Our main aim was to investigate the properties of critical point 
of the dynamical system which play crucially important roles for this interacting model. We have found 
only one critical point which is physically justified. Next we consider small fluctuation about the critical 
point {ucriti Vcrit) and found the eigen values of the corresponding Jacobian matrix for a = 1. It has been 
observed that the all eigen values are negative for all physical choices of parameters and a stable attractor 
scaling solution is obtained. 

In figure 1, we plot dimensionless parameters u and v with respect to x = In a. We notice that u converges 
toward it = 1 in the region < n < 4 (in particular, for n = 1,2,3), also v decreases simultaneously 
and keeps negative sign in near future. This signifies the dark energy dominating nature of the universe 
which means the energy flows from all dark matter to dark energy. Also if we change the values of a, 
then it may be seen that nature of the dimensionless quantities u and v do not sensitively depend on the 
change of a. Figure 2 depicts that the ratio v/u — w vmcg gets the same nature with the variation of n. 
Since in a spatially flat universe, the physically meaningful range of u cr n is < u cr n < 1 leads to the 
condition < c < 3 ( 1+1 "dm)U+a) " along with n < mm{3(l + Wdm)0- + a )j4}, for existence of the critical 
point. This condition of c represents there is an energy transfer from to dark matter to VMCG. The phase 
diagram in u-v space (figures 3,4) shows the attractor solution for n = 1, 3. In the above range, the critical 
point (u cr it, Vcrit) is stable. Thus the present state and the future evolution do not depend sensitively on 
the choice of the initial condition. Our model cannot cross phantom divide, which is consistent with the 
previous work [26], since we have chosen power law form of B(a). If we choose another form of B(a), then 
the VMCG model may cross the phantom barrier, which is more difficult than this one. We may consider 
this type of problem in near future. 

Moreover the expansion of the universe is governed by a power law form around the critical point. Hence 
the expansion will go on forever with an ever increasing rate. We have also studied the evolution of the 
deceleration parameter q(x) = — ^fr = — 1 + | (1 + v(x) + Wd m (l — u(x) j) for the interacting VMCG model 
with dark matter. The result is shown in figure 5 with different parameters value ranging from n — 1 to 
n = 3. Therefore, one can observed that for n = 1,2,3 the values of q decrease from positive value to 
negative value but cannot reached —1. On the other way, figure 6 shows that q decreases from positive value 
to —1 for n = 0. So VMCG (n ^ 0) generates quiessence scenario for late stage but only MCG (n = 0) 
generates the ACDM model for late stage of the Universe. Also we have obtained the statefinder parameters 
at the critical point. Around the critical point, the natures of the parameters are drawn in figures 7 and 8. 
In all the above figures, we have chosen c = 0.01, A = 0.3, Wd m — 0.01 and a — 0.8. 

Distance measurements of the Universe around the critical point have been discussed. The look back 
time, proper distance, luminosity distance, angular diameter distance, comoving volume, distance modulus 
and probability of intersecting objects before and after the present age of the universe have been calculated 
in terms of z, X and Hq and their progressions are shown in figures 9 -16. For all figures, we have 
assumed Hq = 72km/ s/Mpc. The maximum angular diameter distance has been found for a particular 
value of redshift z. Comoving volume element dVc/dz before and after the present age of the universe 
are respectively drawn in figures 11 and 12. We see that when z > 0, the volume element decreases as z 
decreases and when z < 0, the volume clement increases as z decreases. Distance modulus Dm as a function 
of redshift before and after the present age of the universe have been shown in figures 13 and 14. We see 
that for z > 0, Dm decreases as z decreases but for z < 0, Dm increases as z decreases upto a certain stage 
(about z ~ —0.6) and after that Dm decreases as z decreases. In figures 15 and 16 we draw intersection 
probability P as a function of redshift before and after the present age of the universe respectively. We see 
that for z > 0, P decreases as z decreases and for z < 0, P increases as z decreases. So these are the main 
consequences of the power law form of scale factor. 
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